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Abstract. Recent theoretical progress in calculations of the spectrum 
and decays of charmonium is reviewed. Traditionally, our understand- 
OO ■ ing of charmonium was based on potential models. This is now being 

^ ! replaced by first principles. 

o 
in 

-g^: INTRODUCTION 

The charmonium system is at present among the theoretically best understood 
' hadronic systems. The charm quark mass is large compared to the typical QCD 

scale, Aqcd- The cc bound states are therefore governed by non-relativistic dy- 
^ '. namics. Historically, while the QCD potential was not known from first prin- 

^ I ciples, relatively simple guesses for phenomenological potentials had proven 

quite successful in describing the experimentally measured bound state spec- 
trum of charmonium [Q]. Likewise, for annihilation decays of charmonium, an 
intuitive factorization ansatz, involving the wave function at the origin (calcu- 
lable in potential models) worked quite well, for the most part. 

This model-based theoretical understanding is now being replaced by first 
principles. For charmonium spectroscopy, the progress comes from using lattice 
QCD. Charmonium decays can be treated rigorously with the help of non- 
relativistic QCD (or NRQCD). 

Lattice field theory offers a systematic first principles approach to solving 
QCD. It has been argued by Lepage that charmonium is one of the easiest 
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systems to study with lattice QCD, with the potential of complete control 
over all systematic errors. Finite-volume errors are much easier to control for 
quarkonia than for light hadrons. Lattice-spacing errors, on the other hand, can 
be larger for quarkonia and need to be considered. An alternative to reducing 
the lattice spacing in order to control this systematic error is improving the 
action (and operators). For quarkonia, the size of lattice-spacing errors in a 
numerical simulation can be anticipated by calculating expectation values of 
the corresponding operators using potential model wave functions. They are 
therefore ideal systems to test and establish improvement techniques. Most of 
the work of phenomenological relevance is done in what is generally referred 
to as the "quenched" (and sometimes as the "valence" ) approximation. In this 
approximation gluons are not allowed to spht into quark - anti-quark pairs 
(sea quarks). In the case of charmonium, potential model phenomenology can 
be used to estimate this systematic error. Control over systematic errors in 
turn allows the extraction of Standard Model parameters from the quarkonia 
spectra. 

Non-relativistic systems, like charmonium, are best described with an ef- 
fective field theory, non-relativistic QCD (or NRQCD). It was shown in Ref. 
that the application of NRQCD to the problem of annihilation decays of 
quarkonium leads to a general factorization formula. It reproduces the ear- 
lier {ad hoc) factorization ansatz for S-wave decays while putting it on a firm 
theoretical footing. In the case of P-wave decays, the previous ansatz is modi- 
fied. 

Lattice QCD and NRQCD are introduced in the following subsections. The 
remainder of the talk is organized in two parts. The first part reviews recent 
progress in calculations of the charmonium spectrum based on lattice QCD. 
The second part reviews progress in understanding charmonium decays based 
on NRQCD. 

An Introduction to Lattice QCD 

Lattice Field theory is formulated using the Feynman path integral in Eu- 
clidean space. The quantities that are actually calculated are expectation val- 
ues of Greens functions (Q), which are products of gauge and fermion fields. 
The physical quantities of interest, hadron masses, matrix elements, etc., are 
then extracted from these Greens functions. 

The discretization of space-time (with lattice spacing a) regulates the path 
integral at short distances or in the ultraviolet. A finite volume (of length 
L) is necessary for numerical techniques and also introduces an infrared cut- 
off or momentum-space discretization. The vacuum expectation of a Greens 
function, Q, is defined as: 

{g)= \im lim {g)L,a , {G)L,a = Zl\jvi;Vi,VUGe-''^^ . (1) 
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Zi^^a normalizes the expectation value. I have omitted spin and color indices 
for compactness. The gauge degrees of freedom are written as (path ordered) 
exponentials of the gauge field, A^: 

which makes it easy to maintain gauge invariance. The link fields, U, are SU (3) 
matrices. The (Euclidean) QCD action, 

S = Sg + Sf , Sg = ^j d^xF^^F^' , Sf = j d^x^{x){p+ m)^(x) . 

^ (3) 

is discretized, such that Eq. (|^) is recovered in the the continuum (a 0) 
limit: 

^iat = ^ + 0(a") , n>l . (4) 

I will not go into the explicit formulations of Siat here, but instead refer the 
reader to pedagogical introductions . The most common form for the gauge 
action is Wilson's written in terms of plaquettes - products of U fields 
around the smallest closed loop on a lattice. Wilson's gauge action has dis- 
cretization errors of O(a^). 

For fermions the situation is more complicated. The discretization of 

M = p+m, (5) 

is a sparse, finite dimensional matrix. Two different approaches are in use. In 
Wilson's formulation chiral symmetry is explicitly broken, but restored in 
the continuum limit. The pay-off is a solution of the so-called fermion doubling 
problem. Staggered fermions [|^ keep a f/(l) chiral symmetry at the expense 
of dealing with 4 degenerate flavors of fermions. 

Eq. (|I]) emphasizes that QCD is a limit of lattice QCD. However, in nu- 
merical calculations these limits cannot be taken explicitly, only by extrapola- 
tion. This is feasible, because theoretical guidance for both limits is available. 
The zero-lattice-spacing limit is guided by asymptotic freedom, since the lat- 
tice spacing is related to the gauge coupling by the renormalization group. 
Quantum field theories in large but finite volumes have also been analyzed 
theoretically p. 

In a numerical calculation the limits are taken by considering a series of 
lattices, as illustrated in Figure |l|. While keeping the physical volume (or L) 
fixed, the lattice spacing is successively reduced; then, keeping the lattice 
spacing fixed the volume is increased. The calculation is in the continuum 
(infinite volume) limit once the hadron spectrum or matrix elements of interest 
become independent of the lattice spacing (volume). 
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Figure 1: Illustration of the continuum and infinite- volume limits. 

In practice, however, limitations in computational resources do not permit 
the ideal lattice QCD calculation just described. In particular, the compu- 
tational cost of reducing the lattice spacing naively scales like (L/a)^. (The 
computational cost is really higher, because of numerical problems at smaller 
lattice spacings.) Eq. (|^) illustrates an alternative. By improving the discretiza- 
tion errors in the lattice action (and operators), the continuum limit can be 
reached at coarser lattice spacings than before. Simulations with improved 
actions can come at only a slightly higher computational price. The ideas un- 
derlying improvement were developed some time ago |TU|, |TT|, and have since 
been revitalized |13|, |14H|]. 

If the quark mass is large compared to the typical QCD scale, Aqcd, 
effective theories (such as NRQCD) are most adequate in describing the physics 
11^ . In that case, the lattice spacing cannot be taken to zero. Lattice-spacing 
errors can, however, be systematically reduced by improvement |T7| . 

The problem is now (more or less) set up. I again refer the reader to 
the literature |^ for more details on the organization of typical lattice QCD 
calculations. 



An Introduction to NRQCD (Non-Relativistic QCD) 

In non-relativistic systems, the velocity v of the heavy quark inside the bound 
state is a small parameter; in charmonium f ^ ~ 0.3. The important momentum 
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scales with regard to the structure of the bound state are the heavy quark 
momentum (mv) and its kinetic energy (mw^). Momenta of the order of the 
heavy quark mass (m), or above, are relatively unimportant in the bound state 
dynamics. If v is small enough, then all the different momentum scales are well 
separated, in particular: 

Aqcd , mv'^ , mv <^m (6) 

The most adequate description of such systems is in terms of an effective 
field theory, non-relativistic QCD (or NRQCD). The theory has a cut-off A ~ 
m. The effective lagrangian can be written as an expansion in powers of 1/m. 



The NRQCD lagrangian is |]T6 



-^NRQCD — -^^light + heavy + SC . (7) 

flight is the fully relativistic lagrangian for the light quarks and gluons. £heavy 
is the heavy quark (and anti-quark) lagrangian, 

>Cheavy = (^^Do + ^ + (^tD^ " X , (8) 

where ip and x are 2-component Pauli spinors describing quark and anti-quark 
degrees of freedom. 

Relativistic effects of full QCD introduce corrections to Eq. (P) which ap- 
pear in 6C as local, non-renormalizable interactions, with coefficients that are 
calculable in perturbation theory |]I6|, |T^. In principle, infinitely many terms 
must be considered to reproduce full QCD. In practice, however, only a finite 
number is needed, since every operator scales with a certain power of f , as 
shown in Ref. |T^. These power counting rules (e.g., ip ~ (mv)^^^'^\ D ~ mv, 
etc.) effectively order the terms in the NRQCD lagrangian by powers of v. 

The annihilation of a QQ pair occurs at momenta of order m. This short 
distance physics cannot be treated directly in NRQCD. However, the annihi- 
lation contribution to low-energy QQ —>■ QQ scattering can be incorporated 
in NRQCD by adding local 4-fermion operators to 5£ Q, 

(5/^4_fermion — ^ ^, ^ ^ C^i • (9) 

I 

For example, 

(^i05o)=V^^xxV • (10) 

The fi are again calculable in perturbation theory as expansions in as{m). 



5 



CHARMONIUM SPECTROSCOPY 



Two different formulations for fermions liave been used in lattice calculations 
of these spectra. Lepage and collaborators H, [l^ have adapted the NRQCD 
formalism to the lattice regulator. Several groups have performed numerical 
calculations of quarkonia in this approach. In Ref. [jl8| the NRQCD action is 
used to calculate the charmonium spectrum, including terms of 0{mv'^) and 
O(a^). In addition, this group has calculated the bb spectrum in the quenched 
approximation {uf = 0) |T^ and also using gauge configurations that include 
2 flavors of sea quarks 

The Fermilab group [jr3| developed a generalization of previous approaches, 
which encompasses the non-relativistic limit for heavy quarks as well as Wil- 
son's relativistic action for light quarks. Lattice-spacing artifacts are analyzed 



for quarks with arbitrary mass. Ref. [22| uses this approach to calculate the 



cc (and bb) spectra in the quenched approximation. We considered the effect 
of reducing lattice-spacing errors from 0{a) to O(a^). 

The two groups mentioned above use gauge configurations generated with 
the Wilson action, leaving O(a^) lattice-spacing errors in the results. The 
lattice spacings, in this case, are in the range a ~ 0.05 — 0.2 fm. Ref. [p!5| 



uses 



an improved gauge action (to (9(a^)) together with a non-relativistic quark 
action improved to the same order (but without spin-dependent terms) on 
coarse (a ~ 0.4 — 0.24 fm) lattices. 

The first step in any lattice QCD calculation is the determination of the 
two free free parameters of the theory, the gauge coupling and quark mass, 
from experiment. This is discussed in the following subsection. The results for 
the charmonium spectrum from all groups are summarized in Figure ^ 

The agreement between the experimentally-observed spectrum and lattice 
QCD calculations is respectable. As indicated in the preceding paragraphs, 
the lattice artifacts are different for all groups. Figure ^ therefore emphasizes 
the level of control over systematic errors. I should also note, however, that 
the theoretical errors (from the Monte Carlo integration alone) are still much 
larger than the experimental ones. They also grow with every level of excitation 
considered. 

The first quarkonium results with 2 flavors of degenerate sea quarks have 
appeared 2^, 24 1 with lattice-spacing and finite- volume errors similar to the 
quenched calculations, significantly reducing this systematic error. In Refs. |23, 



|2^ the IP-IS splitting in charmonium is calculated, while Ref. considers 
the bb spectrum. 

Several systematic effects associated with the inclusion of sea quarks must 
still be studied. They include the dependence on the quarkonium spectrum of 
the number of flavors of sea quarks and the sea-quark action (staggered vs. 
Wilson) . The inclusion of sea quarks with realistic light-quark masses is very 
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Figure 2: A comparison of lattice QCD results for the cc spectrum using the 
quenched approximation (n/ = 0). The error bars are statistical only. — : 



Experiment; □: FNAL p|; o: NRQCD !§; O: ADHLM |15 



difficult. However, quarkonia are expected to depend only very mildly on the 
masses of the light quarks. This systematic error has not been included yet 
and should be checked numerically. 

The first (and second) generation of lattice QCD calculations, as described 
here, is focused on the simplest physical quantities, like the low lying states or 
simple (decay) matrix elements, to establish the method. Once first principles 
calculations have been achieved, this technology can and will (given sufficient 
motivation) be used to look at more complicated problems, like higher excited 
states, mixing with glueballs, hybrids, etc. 

Standard Model Parameters from Charmonium 

The first step, the determination of the lattice gauge coupling and quark mass, 
follows from comparing appropriate quantities calculated on the lattice with 
the corresponding experimental measurements. The lattice parameters can 
then be converted to their counterparts in continuum QCD with perturbation 
theory. 

Precise determinations of Standard Model parameters are an interesting 
by-product of lattice QCD calculations of charmonium (and bottomonium) . 
However, the theoretical uncertainties must be reduced by an order of magni- 
tude before they become comparable to the present experimental errors. After 
discussing the determination of the lattice spacing (which sets the scale in 
Figure |^), I will summarize the results for the strong coupling and the charm 
quark mass from charmonium. 
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Table 1: Spin-averaged splittings in the J /if) and T systems in comparison. 
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Determination of the Lattice Spacing, a 

The input gauge couphng sets the lattice spacing, a, which is determined in 
physical units by comparing a suitable quantity on the lattice with its ex- 
perimental value. For this purpose, one should identify quantities that are 
insensitive to lattice errors. In quarkonia, spin-averaged splittings are good 
candidates. The experimentally observed IP-IS and 2S-1S splittings depend 
only mildly on the quark mass (for masses between and mc), as shown in 
table |l[ Figure ^ shows the observed mass dependence of the IP-IS splitting 
in a lattice QCD calculation. The comparison between results from different 
lattice actions illustrates that higher-order lattice-spacing errors for these split- 
tings are small pp], In contrast. Figure ^ shows the hyperfine splitting as 
an example of a quantity that strongly depends on both the mass and the 
lattice action. It would therefore be a poor choice for a determination of the 
lattice spacing. 
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Figure 3: The IP-IS splitting as a function of the IS mass (statistical errors 
only) from Ref. □: 0{a^) errors; x: 0{a) errors. 
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Figure 4: The hyperfine splitting as a function of the IS mass (statistical errors 



only) from Ref. □: 0{a^) errors; x: 0{a) errors. 



The Strong Coupling, as 

Within the framework of lattice QCD the conversion from the bare to a renor- 
malized coupling can, in principle, be made non-perturbatively An al- 



ternative is to define a renormalized coupling through short distance lattice 
quantities . The size of higher-order corrections associated with the above 
defined coupling constant can be tested by comparing perturbative predictions 
for short-distance lattice quantities with non-perturbative results . 

At this point the relation to the MS coupling is known to 1-loop, leading 
to a 5 % uncertainty. It has recently been calculated to 2-loops in the 
quenched approximation (no sea quarks, Uf = 0). The extension to nj ^ 
will significantly reduce the uncertainty due to the use of perturbation theory. 

Sea Quark Effects. Calculations that properly include all sea-quark effects 
do not yet exist. If we want to make contact with the "real world" , these effects 
have to be estimated phenomenologically or extrapolated away. 

The phenomenological correction necessary to account for the sea-quark 
effects omitted in calculations of quarkonia that use the quenched approxima- 
tion gives rise to the dominant systematic error in these calculations P5| , 
Similar ideas were used to correct for sea-quark effects in early calculations 
of quarkonia spectra from the heavy-quark potential calculated in quenched 
lattice QCD g. 



By demanding that, say, the spin-averaged IP-IS splitting calculated on 
the lattice reproduce the experimentally observed one (which sets the lattice 
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spacing, a~^, in physical units), the effective coupling of the quenched poten- 
tial is in effect matched to the coupling of the effective 3 flavor potential at 
the typical momentum scale of the quarkonium states in question. The dif- 
ference in the evolution of the zero flavor and 3,4 flavor couplings from the 
effective low-energy scale to the ultraviolet cut-off, where is determined, is 
the perturbative estimate of the correction. 

For comparison with other determinations of a^, the MS coupling can be 
evolved to the Z mass scale. An average [^] of Refs. ||2^, ^ yields for ct^ from 



calculations in the quenched approximation: 

ag^(m2) =0.110 ±0.006 . (11) 

The experimental error in this determination from the quarkonium mass split- 
ting is much smaller than the theoretical uncertainty and does not contribute 
to the total. 

The phenomenological correction described in the previous paragraphs has 
been tested from first principles in Refs. [^, |2^, |4|. All groups calculate 
quarkonium splittings with and 2 flavors of sea quarks. After extrapolating 
to the physical 3 flavor case and evolving the coupling to mz, Refs. [|3], |2^ 
find for the strong coupling from charmonium 

ag^(mz) = 0.111 ±0.005 (12) 

in good agreement with the previous result in Eq. (|11|). The total error is now 
dominated by the rather large statistical errors and the perturbative uncer- 
tainty. 

At present, the result of Ref. ||2^ has the smallest statistical and systematic 
errors for the strong coupling (in this case from the bb spectrum): 

ag^(mz) =0.115 ±0.002 . (13) 

Phenomenological corrections are a necessary evil that enter most coupling 
constant determinations. In contrast, lattice QCD calculations with complete 
control over systematic errors will yield truly first-principles determinations of 
as from the experimentally observed hadron spectrum. 

At present, determinations of from the experimentally measured quarko- 
nia spectra using lattice QCD are comparable in reliability and accuracy to 
other determinations based on perturbative QCD from high energy experi- 
ments. They are therefore part of the 1994 world average for as [^. The 



phenomenological corrections for the most important sources of systematic er- 
rors in lattice QCD calculations of quarkonia are now being replaced by first 
principles, which will significantly increase the accuracy of as determinations 
from quarkonia. In particular, the systematic errors associated with the inclu- 
sion of sea quarks into the simulation have to be checked. 
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The Charm Quark Mass 



Because of confinement, the quark masses cannot be measured directly, but 
have to be inferred from experimental measurements of hadron masses, and 
depend on the calculational scheme employed. In lattice QCD quark masses are 
determined non-perturbatively, by tuning the input lattice quark mass (wig*) 
so that, for example, the experimentally observed J/ip mass is reproduced by 
the calculation. 

Phenomenologically useful quark masses are the perturbatively defined pole 
and MS masses, which the bare lattice mass can be related to by perturbation 
theory: 

^Poie^^poie^gt mf{mQ) = Z^m'^^ . (14) 

Of course, as always, all systematic errors arising from the lattice QCD cal- 
culation need to be under control for a phenomenologically interesting result; 
in particular, the systematic error introduced by the (partial) omission of sea 
quarks has to be removed. The short-distance corrections that introduced the 
dominant uncertainty to the as determination from quarkonia are absent for 
the pole mass determination, because this effective mass does not run for mo- 
menta below its mass. An analysis of the 6-quark mass from the hh spectrum 
with and without sea quarks is consistent with this estimate [pif . 



Ref. [0 analyzes the the charm quark mass from the charmonium spec- 



trum with the preliminary result, m^°'*^ = 1.5(2) GeV. 

The MS mass for the charm quark has also been determined from a com- 
pilation of D meson calculations in the quenched approximation with 
m^^(2GeV) = 1.47(28) GeV. The error includes statistical errors from the 
original calculations and the perturbative error. However sea-quark effects 
cannot, in this case, be estimated phenomenologically, leaving this system- 
atic error uncontrolled. 



CHARMONIUM DECAYS 

Historically p4| , charmonium annihilation decays were treated with a factor- 
ization ansatz, which divides the decay rate into a short distance, perturbative 
part, and a long distance, non-perturbative part, usually parametrized as the 
wave function at the origin, R{0). This ansatz worked quite well for S-wave 
decays, even after including radiative corrections at next-to-leading order p5| | . 
For P-wave decays, the long distance piece was identified as the derivative of 
the wave function at the origin, R'{0). However, the radiative corrections were 
found to have infra-red divergences |^6|] for J = 1 P-wave states already at 
leading order, and for J = 0,2 P-wave states at next-to- leading order, pre- 
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venting reliable theoretical predictions. This ansatz also did not allow for a 
systematic inclusion of higher Fock states, or relativistic corrections. 

If the effective field theory framework of NRQCD is used to describe anni- 
hilation decays, it was shown that a general factorization formula holds for 
the decay rates. The annihilation decay rate of a quarkonium state A can be 
written as 

r(A) = 2 Im(A|5/:4-fcrmion|^) • (15) 

Using Eq. @ this gives 



nA) = T.;;^,{m\A) . (16) 

i 

The coefficients Fi are the imaginary parts of the in Eq. (M). Because of the 



power counting rules of Ref . |Tj] , the matrix elements scale with some power 
of the velocity, 

{A\0,\A) r-. , (17) 

which usually increases as higher dimensional operators are considered in the 
sum of Eq. (|^). Effectively, this approach expresses the quarkonium decay 
rates as expansions in the the short distance parameter as{m) and the long 
distance parameter v"^. The desired accuracy of the theoretical prediction thus 
serves as a truncation criterion for the (infinite) sum in Eq. (|T6|). 

The matrix elements in Eq. ([T6|) are calculable from first principles using 
lattice NRQCD. These calculations are in progress by the ANL group [^. In 
the meantime, the matrix elements can also be extracted from experimental 
measurements of charmonium decay rates. 

It was shown in Ref. 0] that heavy-quark spin symmetry and the vacuum- 
saturation approximation reduces the number of independent matrix elements 
that have to be determined non-perturbatively or phenomenologically. 

The matrix elements of the 4-fermion operators, which contribute to the 
charmonium decays into light hadrons can be simplified using the vacuum sat- 
uration approximation. This approximation is valid in NRQCD up to relative 
order f ^. For example, 

M^xx^i^lvc) = Ko|xVl^c)P(i + o{v')) . (18) 

This relates the matrix elements appearing in electro-magnetic charmonium 
decays to those in strong decays (into light hadrons). 

Heavy-quark spin symmetry relates matrix elements of states in the same 
radial and orbital levels but different spins. For example, 

e* ■ (OlxVV'l J/^) = (0|xVl^c)(l + 0{v^)) . (19) 

The following two subsections discuss the application of this approach to 
S- and P-wave decays, using specific examples. 
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S-Wave Decays 

I discuss the theoretical knowledge of S-wave decays using r]c — > 77 as an 
example. At leading order in f ^ the rate is 

nVc~^ll) = ^^^^\Rs? , (20) 

where by spin symmetry (see Eq. (|19|)) Rs can be taken as i?^^ or the spin- 
averaged combination of Rr,^ and i?^, with 

^=\/^(0|xVl^c) . (21) 

The other rjc and J/ip decays differ from Eq. (^) in the short distance coef- 
ficient but, using Eqs. (|1^) and (p!9|), have the same long distance parameter, 
Rs- Identifying Rs with the wave function at the origin, -Rs(O), leads back to 
the old factorization ansatz. The short distance coefficients, Fj, are known to 
next-to-leading order in a^- 

At next to leading order in f^, the decay rate becomes 

1 c \ n ' 



F^-v( So) ,— — ,o G'wf Si 



- 77) = \Rr,f + ^^^^ ReiRsV'Rs) . (22) 

Now the difference between i?^^ and Rs, which is of order (see Eq. (|l^)) 
has to be taken into account for consistency. V^-R^ can be interpreted as the 
laplacian of the wave function at the origin. 



{0\xKD')^|J\r^^) = V'Rs (1 + 0{v')) . (23) 

The short distance coefficients, Gi, are known to leading order in as- 

The long distance parameters, Rrj^, R^ and 'V^Rs can be calculated 

in lattice NRQCD, estimated using potential models or extracted from a 

phenomenological analysis of experimental data. The relative errors (from 

the truncation of the perturbative and non-relativistic series) are ^^(mc)^, 

as{mc)v^ and v"^, adding up to <15 %. 

It is not inconceivable that at least for the electro-magnetic S-wave decays 

some (or all) of the higher order calculations will be performed, reducing the 

theoretical error accordingly. 

P-Wave Decays 

The structure of P-wave decay rates is more complicated than that of S-waves. 
I shall start with a simple example first, the decay Xcj 77- 

nxcj - 77) = ^^^^ , (24) 



13 



where again, by spin symmetry R'p can be taken from the matrix elements of 
any of the P-wave states (or their spin-averaged combination). For example, 

1 / 27 

{Q\x\-^D . <t)Mxco) = (1 + Oiv')) . (25) 

The obvious interpretation of R'p as the derivative of the wave function at the 
origin, R'p{0), again connects this formalism to the old factorization ansatz. 
The coefficients, F^^, are known to next-to-leading order in a^. The relative 
error is v"^ or ~ 25 % of the decay rate. 

The situation is different for strong decays of P-wave states. Taking as an 
example he —>■ LH (light hadrons), Eq. (|1^) gives at leading order in f^, 

p np \ p (I c \ 

T{h, LH) = ^V^(Oi) + ^^(Os) , (26) 



where 



(Oi) = {K\i;\]-D ■ (t)xxKId • cT),lj\h,) (27) 



2 "^'^ '2 
R'pf (l + Oiv^)) 



271 



and 



{Os) = {Kl^/j^T^x^T^^lhc) ■ (28) 

The two matrix elements (Oi) and (Og) enter at the same order in w^, because 
they are both suppressed by f ^ with respect to the leading order S-wave matrix 
elements. Oi is a dimension eight operator; the two powers of D give the 
suppression. Cg is of dimension six, but its matrix element picks the \ccg) Fock 
state, where the cc pair is in a color-octet ^5*0 state. The dominant Fock state 
of charmonium is the color-singlet | cc) . 

\A) = ^,,\cc) + '$,sg\ccg) + ... , (29) 

with \E'ccs ~ C'(f ). This gives the v"^ suppression. 

This departure from the old factorization ansatz solves the problem of 
infrared divergences, and thus leads to a consistent treatment of this case 
similar to the other charmonium decays. 

The short distance coefficients, Fi, are known to next-to- leading order, 
while the Fg's are only known to leading order in a^. At present, the dominant 
relative error is thus as and v^, adding up to ~ 40 %. 

It should be possible to extend the (perturbative) calculations of P-wave 
decays in this frame work to the same order as what is presently available for 
S-wave decays, reducing the theoretical error to <15 % (assuming knowledge 
of the relevant matrix elements). 

A comparison of theory and experiment for P-wave decays shows fair agree- 
ment, albeit with rather large errors from both sides p8| . 
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CONCLUSIONS 



Quarkonia were, upon their discovery, called the hydrogen atoms of parti- 
cle physics. Their non-relativistic nature justified the use of potential models, 
which gave a nice, phenomenological understanding of these systems. This 
phenomenology is at present useful to control systematic errors in lattice 
QCD calculations of the charmonium spectrum. However, we are quickly mov- 
ing towards truly first-principles calculations of quarkonia using lattice QCD, 
thereby testing QCD non-perturbatively. In this sense, quarkonia are still the 
hydrogen atoms of particle physics. Precise determinations of the Standard 
Model parameters, Ofg, rUc (and m^), are by-products of this work. 

Still lacking for a first-principles result is the proper inclusion of sea quarks. 
The most difficult problem in this context is the inclusion of sea quarks with 
physical light quark masses. At present, this can only be achieved by extrap- 
olation (from rriq ~ 0.3 — 0.5ms to rriu^d)- If the light quark mass dependence 
of the quarkonia spectra is mild, as anticipated, the associated systematic er- 
ror can be controlled. First-principles calculations of quarkonia could then be 
performed with currently available computational resources. 

The present theoretical status of charmonium annihilation decays is rather 
promising. The frame-work developed in Ref. |^ leads to a systematic expan- 
sion in as{mc) and v'^, with controllable uncertainties. Until first-principles 
calculations of the non-perturbative matrix elements become available, this 
formalism can still be tested phenomenologically, using experimental data. 
Theoretical predictions for the decay rates should be available with uncertain- 
ties of <10 %, in most cases before the Tau/Charm factory turns on. 

It is conceivable, that by the time a Tau/Charm factory turns on, the the- 
ory of charmonium will be solidly based upon first principles with accurate 
predictions for the spectrum and decays of the low-lying states. We will have 
moved to the next stage in theoretical (first-principles) calculations concern- 
ing, for example, the properties of hybrid states or mixing with glueballs. Ex- 
perimental information gathered at the Tau/Charm factory and earlier experi- 
ments will then give us precision tests of perturbative and non-perturbative 
QCD in the charmonium system. 
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